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Abstract 

In these lectures we first briefly review the main observational facts which 
imply that most part of matter in the Universe is not visible and some recent 
intriguing experimental data which would point to a significant contribution 
to J7 due to a cosmological constant. We subsequently discuss some particle 
candidates for dark matter, with particular emphasis for the neutralino. We 
present the main properties of this particle, also in the light of the most recent 
experimental results in direct search for relic particles; furthermore, we discuss 
the perspectives for their indirect searches. 



1. Dark Matter and Cosmology 

The two pillars of standard cosmology are the Einstein equations 

Rm. - ^T^St^u = SttGT^, + Ag^, , (1) 

(where R^u is the curvature (Ricci) tensor, 71 is its trace, T^,^ is the stress-energy 
tensor, A the cosmological constant) and the Robertson- Walker metric 
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ds'^ = dt^ - YCit) — + r'^de'^ + sin^ Odcj)'^] , (2) 

\l — kr'^ I 

where R = R(t) is the cosmic scale factor and t, r, 6, 6 are the comovine coordinates 

We recall that Eq. follows from the assumption that the distribution of matter 
and radiation in the Universe is isotropic and homogeneous. When the coordinates 
are appropriately rescaled, the values k = +1, 0, —1 define the space curvature to be 
positive, zero and negative, respectively. Useful standard quantities for the description 
of the expanding Universe are the Hubble parameter (expansion rate of the Universe) 
H = R/ R and the deceleration parameter q = — R ■ R/R^, whose values at present 
epoch are denoted as Hq and qo. 

By combining Eqs. (|l]) - one obtains the Friedmann equation 

nm + ^A- k/iR^R") = 1 , (3) 

where flm is the ratio of the average matter-energy density p to the critical density 
Pc = 3H^/(87rG'), Qm = p/Pc, and = A/(3H^). The critical density is given by 
Pc = 1.9 ■ 10^^^ h"^ g cm~'^, when the Hubble constant is parametrised as follows: 
h = Ho/(100 km s'^ Mpc^^). 

From Eqs. (|T]) - it also follows that at present time, when the radiation 
contribution to Q can be set to zero, the value of the deceleration parameter is given 
by 

go = - • (4) 

The features of the evolution of our Universe depend on the actual values to be 
assigned to the cosmological parameters previously defined. In what follows we briefly 
summarize some of the main observational data about these parameters. 



2. Age of the Universe and Expansion Rate 

The evaluation of the present age of the Universe, to, depends on the expansion 
rate and on the various components of Q; therefore measurements of to and of the 
Hubble constant provide information on the size of Qm and (see, for instance, 
Refs. 11). 

A recent determination of to provides the value: to = 11-5 ± 1.3 Gyr Q, with a 95% 
CL. lower bound of 9.5 Gyr. Recent measurements of the Hubble constant by the 
Hubble Space Telescope Key project give the following range i'i: Hq = 73 ± 6 (stat) ± 
8{syst) km s~^ Mpc~^. In view of the still persisting sizeable spread in the h values, 
due to a host of independent measurements (for a review, see for instance i) in the 
following we will, conservatively, consider a rather wide range: 0.5 < h < 0.9. 
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If, for sake of illustration, we take ~ 0.7, it is easy to show that to ~ 11-5 Gyr 
would require either < 0.3 or ~ 0.5, according to whether or not we allow for 
a non- vanishing cosmological constant, such that Q\ = 1 — Qm- 

3. Observational Evidence for Dark Matter 

The parameter Q^n niay be derived from astronomical determinations of the aver- 
age mass-to-light ratio M/L for various astrophysical objects (see, for instance, Ref. 
B for an updated review) 

_M L 

L Pc 

where L is the luminosity density: L = 1.6-10® h Lq Mpc~^. Using pc = 2.8-10^^ h'^ Mq 

Mpc~^, one gets 

_ h-' M/L 
" 1500 Me/Lg ■ 

From the M/L ratio measured in stars, M/L ~ (3 — 9)Mq/L0 we obtain 

0.002 h-^ < n^,s < 0.006 h-^ . (7) 



3.L Rotational curves of spiral galaxies 

Presence of dark matter in single galaxies is apparent from the flatness of the 
plot of the rotational velocities versus the galactocentric radius, well beyond the 
distribution of the visible matter. From these measurements one derives M/L ^ 
70 M0/L0(_R/iaio/lOO kpc). Then we obtain 

Qhalo-^ 0.05 hRhaio/lOOkpc, (8) 

a result which, compared to the range (0) for fl^is, is indicative of the presence of 
dark matter at the level of single galaxies. 

3.2. Clusters of galaxies 

Existence of dark matter at the level of clusters of galaxies may be established by a 
number of methods: X-ray emission by hot gas in intra cluster plasma, measurements 
of velocity dispersion and strong gravitational lensing. 

We report here some results derived from measurements of the X-ray emission liS. 
For rich clusters of galaxies one finds M/L = (300 ± 100) HMq/Lq, which gives 

^^dn^ter^ 0.2 ±0.07. (9) 
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The baryonic content fib is established to be: ~ 6 /i '^/^% for gas and > 4% for stars; 
then 



fib 



fi 



> 0.04 + 0.06/1"^/^ (10) 



^cluster 

On the other hand, the Big Bang nucleosynthesis provides the following estimate 
for fife l§: 

0.009 < fifo/i' < 0.02 , (11) 

or, taking 0.5 < h < 0.9, 0.01 < fib < 0.08. Combining this result with Eq.(0) we 
obtain 

^cluster < 0.2 - 0.4 . (12) 



3.3. Large-scale Velocity Flows 

Distribution of matter over large scales may be inferred from the peculiar motion of 
the gravitationally-induced inflow. Let us consider the case of the Local Supercluster, 
centered near the Virgo cluster; we stay on the edge of this cluster, at a distance of R 
= 12 h^^ Mpc. The radial inward peculiar velocity averaged over the surface is given 
by@ 

= ^HoR5„fi°^ (13) 

where fi°'^ represents the factor due to expansion and 6m denotes the mass contrast, 
related to the contrast in galaxy counts {Sglhj the relation 6m = 6g/b, b being the 
bias factor. Using the observational values M: v = 200 ± 25 km s~^, Sg = 2.3 ± 0.7, 
and taking 6 ~ 1, one obtains 

fi^^o.ii°:J5. (14) 



3.4- A few first conclusions 

From the previous evaluations for fi„js and fi^, and from the range for fib in 
Eq.([Tl|) we derive the following important points: 



• a large amount of matter in the Universe is not visible; 

• some of this dark matter is baryonic; 

• a significant amount of dark matter is non- baryonic. 
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As usual we divide particle dark matter into two categories: Hot Dark Matter 
(HDM) and Cold Dark Matter (CDM), according to whether the particles are rela- 
tivistic or non-relativistic at their decoupling from the primeval plasma. 

4. Singling out different contributions to Q 

We briefly report now some important observational results and analyses which 
provide further clues toward a determination of the various components to Q. The 
two main issues are: what, if any, is the size of f^A and what is Qm made of? 

4.I. Formation of cosmological structures 

A standard method for testing different dark matter models is to compare the 
power spectrum of the density fluctuations P{k) with measurements of the Cosmic 
Microwave Background Radiation (CMBR) anisotropy and measurements of the two- 
point correlation function in galaxy surveys illfl'i. We recall that P{k) is the Fourier 
transform of the correlation function between the density contrasts at two different 
points in space. It is customary to assume for P{k) a power law, i.e. P{k) oc fc". The 
Harrison-Zel'dovich spectrum corresponds to n = 1. In the past, typically, the best 
performing cosmological models turned out to fall into the following categories: 

• HCDM model characterized hj Q = 1 and the following composition: fib — 0.1, 
fli, ~ 0.2, ficDM — 0.7; h ~ 0.5, where fib, and ficDM are the contribution 
due to baryons, neutrinos (as HDM particles) and to CDM particles, respec- 
tively. 

• ACDM model with fi = 1 and ficoM ^ 0.3, fiA ^ 0.7; /i ~ 0.7 - 0.8. 

• TCDM model (= tilted CDM model) with a power spectrum P{k) oc k^'^ and 
fi = 1 = ficDM = 1; ^ 0.5. 



4.2. Evidence for Qa ^ ? 

Recent measurements of high-redshift supernovae of type la point to an 

important contribution to fi due to A, with a relatively small contribution from 
fim- These data appear to be complementary to those derived from measurements 
of the CMBR t3. The joint use of these two sets of data, together with some other 
observational data, singles out the following ranges: fi^ = 0.24 ± 0.10 and fiA = 
0.62 ± 0.16 0. These data and analyses are of the utmost interest for their potential 
implications, although a number of points need further clarification i. 

One further hint for a rather low value of fim is provided by the time evolution of 
the number density of clusters. In Ref.lll observational data on cluster abundance in 
the redshift range < 2; < 1 is used to derive the estimate fim — 0.2^^0;]^^. 
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Though a cautionary attitude is in order here, it is important to remark that 
all these different ways of determining point to a relative small value for this 
quantity: Q^a < 0.3 - 0.4. This feature, if confirmed by further observational data, 
has profound implications for the phenomenology related to DM particle candidates, 
as will be illustrated in the following. 

5. Particle candidates for dark matter 



5.1. Baryons 

As was already noticed in Sect. Ill, some contribution to DM is provided by 
baryons. This conclusion is drawn from the fact that the Big Bang nucleosynthesis 
provides the estimate 0.009 < Qfyh"^ < 0.02, which, together with Eq. (|^, implies 

^ ^vis- 

Direct search for non-luminous baryonic dark matter, under the form of microlens- 
ing objects, has been undertaken since the seminal paper of Ref.0. Recent results 
are reviewed in The main features of the data may be summarized as fol- 

lows. The EROS Collaboration @ excludes that microlensing objects of masses in the 
range (5 ■ 10^^ — 10^^) M© may make up more than 20% of the halo density, whereas 
the MACHO Collaboration ^ delimits a likelihood contour (at 95% C.L.) for masses 
~ (10~^ — 1) Mq with a best-fit value for the halo fraction around 0.5. For further 
details see, for instance, Refs. 11013. For a critical view of the microlensing events, 
see Ref.il. An interesting scenario related to baryonic dark matter is the one depicted 
in Ref.i. 



5.2. Non-baryonic DM candidates 

Particle physics offers a large variety of particles, which would have decoupled 
from the primeval plasma at the time (freeze-out time), when the interaction rate 
became smaller than the cosmic expansion rate; these particles would then be floating 
around in our Universe as relics. 

These fossil particles would or would not significantly contribute to the average 
cosmic density depending on their actual number density and mass. The most obvious 
example is provided by light fossil neutrinos, whose relic abundance may be easily 
evaluated (see, for instance Ref.i) and turns out to be 

S^^^^^i^, (15) 

where the sum is over the neutrino flavours (for each flavour, neutrino and antineu- 
trino are counted together). Therefore the relevance of these fossil neutrinos for the 
Universe matter density depends on whether their mass is of order of a few eV or 
much smaller (provided neutrinos are massive at all). 
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Possible indications for non-vanishine neutrino masses are from: 1) solar neu- 
trinos atmospheric neutrinos IISS, and 3) the LSND experiment 0. All 
these data refer to oscillation measurements, and then are not sensitive to individ- 
ual neutrino masses, but only to differences in their squares. Typically, atmospheric 
neutrino experiments give Am^ ~ (2 6) x 10~^eV^; solar neutrino experiments 
can be explained either by a Am? ~ 10~^eV^, in case of matter-enhanced oscilla- 
tions, or by a Am? ~ 10^^'' eV^, in case of vacuum oscillations; LSND data suggest 
0.2 eV^ < Am^ < 2eV^. These results may be compatible with a sizeable value of 
relic abundance, ~ 0.2, such as preferred by some calculations of cosmological 
structures. However, if taken at their face values, they only imply > 0.02 (0.5//i)^ 
i. 

A different candidate is the axiom whose motivation in particle physics is related 
to the strong CP-problem in QCD M. A discussion of this candidate is beyond the 
scope of the present report; for a comprehensive review on its possible cosmological 
relevance and on the experimental efforts for detecting it as a relic particle, see for 
instance Ref. 0. 

Among the particle candidates for CDM the most favorite one is certainly the 
Lightest Supersymmetric Particle (LSP), under the condition that it is weakly inter- 
acting. This candidate is discussed in the following section. 



6. Supersymmetric dark matter particles 

In order to be a dark matter candidate a particle has to be weakly interacting 
and stable (or, at least, long lived on cosmological time-scales). A very interesting 
perspective for such a candidate is offered, in the framework of supersymmetric theo- 
ries with R-parity conservation, by the lightest susy particle (provided this is indeed 
weakly interacting). Different candidates have been considered: the neutralino H 
or the sneutrinoij in gravity mediated models, the gravitino0 or messenger field^ 
in gauge mediated theories, the axinoE^, stable non-topological solitons (Q-balls)E3, 
heavy non-thermal relic£ll or others. 

Among the different candidates, the most promising one turns out to be the 
neutralino, since its relic abundance may be sizeable, at the level required to explain 
the CDM content of the Universe and, at the same time, its detection rates may be 
accessible to experimental searches of different kinds. 

The phenomenology of neutralino dark matter has been studied extensively in the 
Minimal Supersymmetric extension of the Standard Model (MSSM) e3. This model 
incorporates the same gauge group as the Standard Model and the supersymmetric 
extension of its particle content. The Higgs sector is slightly modified as compared 
to that of the Standard Model: the MSSM requires two Higgs doublets Hi and H2 
in order to give mass both to down- and up-type quarks and to cancel anomalies. 
After electroweak symmetry breaking, the physical Higgs fields consist of two charged 
particles and three neutral ones: two scalar fields {h and H) and one pseudoscalar (A). 
The Higgs sector is specified at the tree level by two independent parameters: the 
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mass of one of the physical Higgs fields and the ratio of the two vacuum expectation 
values, usually defined as tan/5 = < H2 > / < Hi >. The supersymmetric sector 
of the model introduces some other free parameters: the mass parameters Mi, M2 
and M3 for the supersymmetric partners of gauge fields (gauginos), the Higgs-mixing 
parameter /z and, in general, all the masses of the scalar partners of the fermions 
(sfermions) and all the trilinear couplings which enter in the superpotential. In the 
MSSM it is generally assumed that the gaugino masses are related by expressions 
induced by grand-unification. Specifically, the two parameters which are relevant 
for neutralino phenomenology are linked by: Mi = (5/3) tan^ ^vy^2- The other 
usual assumptions are that all slepton mass parameters are taken as degenerate to a 
common value mo and that all the trilinear couplings are vanishing except the ones of 
the third family which are set to a common value A. In summary, the free parameters 
of the model are six: M2, /i, tan/?, m^, mo and A. 

The neutralinos are four mass-eigenstates defined as linear superpositions of the 
two neutral gauginos (7 and Z) and the two neutral higgsinos {Hi and H2) 

X = ai7 + + a^Hi + a4H2 . (16) 

The lowest-mass eigenstate may play the role of the lightest supersymmetric particle 
in the MSSM, and may then constitute the dark matter candidate in this model. It 
will be called the neutralino tout-court and its mass denoted by m^. To classify 
the nature of the neutralino it is useful to define a parameter P = af + a^; the 
neutralino is called a gaugino, when P > 0.9, a higgsino when P < 0.1 and mixed 
when 0.1 < P < 0.9. 

The low-energy MSSM scheme is a phenomenological approach, whose basic idea 
is to impose as few model-dependent restrictions as possible. At a more fundamental 
level, it is natural to implement this phenomenological scheme within the supergrav- 
ity (SUGRA) framework One attractive feature of the ensuing model is the con- 
nection between soft supersymmetry breaking and electroweak symmetry breaking, 
which would then be induced radiatively. Usually, the low-energy phenomenology of 
SUGRA theories constitutes a subset of the susy configurations which are considered 
in the MSSM. A typical characteristic of SUGRA models is in fact the presence of 
relatively strong correlations among the low-energy parameters, correlation which 
is absent in the MSSM. In this report we will discuss neutralino dark matter phe- 
nomenology in the less constrained MSSM model. Results for SUGRA schemes can 
be found in the papers listed in Ref.@ 

6.1. Neutralino relic abundance 

Neutralinos decouple from the hot plasma in the early Universe when they are no 
longer relativistic. Their present abundance is calculated by solving the Boltzmann 
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equation for the evolution of the density of particle specie^ and can be written as: 



where 



^xh' = C '-^^^ — (17) 



*° - S.7- lO-i^GeV" 



0.264 peMp 

In the previous Eqs., g*{Tf) and g*g(Tf) denote the effective number of degrees 
of freedom for the energy density and for the entropy density, respectively, evaluated 
at the freeze-out temperature T/; (crann^'r)int is the neutralino pair annihilation times 
the pair relative velocity, averaged over the neutralino thermal density distribution, 
integrated from the freeze-out temperature down to the present temperature; Sq 
denotes the present entropy density, pc is the critical density and Mp is the Planck 
mass. 

The critical quantity to be evaluated is the neutralino annihilation cross section, 
which, depending on the neutralino mass, can get contributions from the following 
final states : i) fermion-antifermion pair, ii) pair of neutral Higgs bosons, iii) pair of 
charged Higgs bosons, iv) one Higgs boson-one gauge boson, v) pair of gauge bosons. 
The diagrams contributing to the final state i) are Higgs-exchange Z-exchange dia- 
grams in the s-channel, sfermion-exchange diagrams in the t-channel. For the other 
final states, the contributions come from Higgs and Z-diagrams in the s-channel, 
and either neutralinos or charginos exchange in the t-channel, depending on the elec- 
tric charges of the final particles When the mass of the neutralino is close 
to the mass of another supersymmetric particle, the process of co-annihilation 0@ 
can be present. In this case, the calculation of the annihilation cross section and of 
its thermal average has to take into account a large number competing interactions 
among the neutralino and its close-mass particles. In some special instances the relic 
abundance may be affected by co-annihilation effects by a sizeable amount^. 

In Fig. 1 we show il^h"^ as a function of the neutralino mass 0. We present 
here, as well as in the following, all the results in terms of scatter plots which are 
obtained by varying the supersymmetric parameters inside wide ranges. Naturally, 
the parameter space is constrained by experimental limits on Higgs bosons and su- 
persymmetric particles searches (for recent updates, see ii) and by measurements of 
rare processes, whose theoretical predictions are quite sensitive to supersymmetric 
contributions. At present, the most important is the decay b ^ s + 700. 

The scatter plot of Fig. 1 is obtained by varying the supersymmetric param- 
eters in the following ranges: 20 GeV < Ms < 500 GeV, 20 GeV < < 
500 GeV, 80 GeV < tua < 1000 GeV, 100 GeV < mo < 1000 GeV, -3 < A < 
+3, 1 < tan/3 < 50. 

The figure shows only those configurations which provide a value for the relic 
abundance lower than a cosmological upper bound which has been conservatively 
set as Q^h"^ < 0.7. The susy configurations which entail larger values of ^^h"^ are 
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excluded by the lower limit on the age of the Universe 0. This constraint is rather 
restrictive on the susy parameter space. 



6.2. Detection of relic neutralinos 

Relic neutralinos would act as CDM during the process of galaxy formation. It 
is therefore conceivable that they may constitute all or part of the dark matter halo. 
These neutralinos would be clustered in the Galaxy and hence possess a matter density 
distribution and a velocity distribution which depend on the dynamics of the Galaxy 
formation and evolutionEd. 

Many different models have been discussed in the literature for the matter density 
distribution p{r) (see for instance Refs. HSIl). This field is in rapid expansion, due 
to the hieh resolution simulations now at hand to investigate the structure of single 
galaxies E^. In particular, these studies are expected to pin down the nature of the 
cuspy behaviour which appears to occur near the galactic center. Another very inter- 
esting possibility which is currently being investigated is that the halo could present 
a clumpy distribution of dark matter together with or in alternative to a smooth 
distribution. The uncertainties in the shape profile, combined with experimental un- 
certainties and the possibility that some fraction of the dark halo is made of baryonic 
dark matter in the form of MACHOS, implies that the local value of the matter den- 
sity pi = pItq) is somewhat uncertain. At present, its most conservative range of 
variability can be set as 0; 0.1 GeV cm^^ < < 0.7 GeV cm~^. 

The quantity pi refers to the total dark matter density of the galactic halo, the 
neutralino local density is, in general, a fraction of it, i.e. p^ = C,pi, with ^ < 1. No 
exact way to evaluate the quantity ^ is currently available. A usual procedure is to 
consider the galactic halo as composed entirely of neutralinos if, on the average in the 
Universe, neutralinos are abundant enough to explain all the CDM which is observed 
at least on the galactic scale. This happens when Q^h'^ is larger than a value (f2/i^)min 
derived from astrophysical observations. In this case it is possible to set ^ = 1. If, 
on the contrary, Q-^h"^ < (f2/;,^)min, neutralinos are not able to explain all the CDM, 
even the one which is needed at the galactic scale, and therefore also locally in our 
Galaxy we expect them to give only a fractional contribution C, to pi. In this case, 
the simplest choice is to set: ^ = fl^h"^ /{flh'^)min- The quantity (fi/i^)mm is estimated 
to lie in the range 0.01 < (fi/i2)mm < 0.2. 

The velocity distribution f {v) of dark matter is usually assumed to be a Maxwellian 
distribution (as seen from the galactic rest frame), with a velocity dispersion fj-ms 
which is directly related to the asymptotic flat rotational velocity Voo as: v^- 



^rms 



3/2voo. In our Galaxy, the rotational velocity appears to be already fiat at the local 
position, and therefore we set foo = Vq. The experimental determination of the local 
rotational velocity is: Vq = 230 ± 50 Km s~^ (90 % C.L.)!!. The distribution is also 
truncated bvan escape velocity fesc which falls in the range: fesc = 450 -r- 650 Km s~^ 
(90 % C.L.)E3Ey. Modifications to the standard Maxwell-Boltzmann form have been 
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examined Hill, but the problem of determining the correct form of the distribution of 
the velocities in the halo has no clear and simple solution at present, both theoretically 
and observationally. Also the possibility that the halo could possess a bulk rotation 
has been considered in the literature^. 

Due to the possibility that neutralinos are present in the halo of our Galaxy, it is 
of great interest, both for astrophysics and particle physics, to search for techniques 
capable of detecting these dark halo particles. Several methods have been proposed 
and in the following we will briefly review the ones which, at present, appear to be 
more promising. 



6.3. Direct detection of relic neutralinos 

The most direct possibility to detect dark matter particles is to look for their 
scattering with the nuclei of a low-background detector!!. The interaction of slow 
halo neutralinos of mass > 25 GeV with a detector produces the recoil of a nucleus 
with energy of the order of few to tens keV. The recoil energy can be measured 
by means of various experimental techniques with different nuclear species, like Ge, 
Nal, Xe, CaFs, TeOsS. The relevant quantity to be calculated and compared with 
the experimental measurements is the differential detection rate 



where Nt is the number of the target nuclei per unit of mass, is the local neutralino 
matter density, v and f{v) denote the neutralino velocity and velocity distribution 
function in the Earth frame {v = \v\). The nuclear recoil energy is given by = 
m^gjf^(l — cos9*)/mN, where 6* is the scattering angle in the neutralino-nucleus 
center-of-mass frame, m^r is the nuclear mass and mred is the neutralino-nucleus 
reduced mass. Finally, da/dEji is the neutralino-nucleus differential cross section, 
which has a coherent contribution due to Higgs- and squark-exchange, and a spin 
dependent one which originates from the exchange of the Z boson and squarks. The 
coherent cross section is usually largely dominant over the spin-dependent one. 

Eq. (|T9|) refers to the situation of a monoatomic detector, like the Ge detectors. 
For more general situations, like for instance the case of Nal, the generalization is 
straightforward. From those experimental data on the nuclear recoil spectrum which 
do not provide a signal-to-background discrimination, upper limits to the neutralino- 
nucleus cross section as a function of the neutralino mass may be set by employing 
Eq.([T9|) Ed. In the case of coherent interaction. Fig. 2 shows, as a solid line, the 
present most stringent upper limiti^ on the quantity ^o"scaiar ^s. m^, where £7^^^?^^^'^°"^ 
denotes the scalar elastic cross section of a neutralino off a nucleon. The astrophysical 
parameters are chosen as: pi = 0.3 GeV cm~^, vq = 220 Km s~^, v^sc = 550 Km s~^ 
and {nh^)mm = 0.01. 

In this plot, we also show a scatter plot of the same quantity calculated in the 
The susy configurations have been varied in the ranges quoted in Sect. 
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6.1. We stress that the cosmological bound f^^^^ — ^-'^ been applied E3. 

In the case of direct detection, a typical signature consists in the annual mod- 
ulation of the detection rate0. During the orbital motion of the Earth around the 
Sun, the change of direction of the rehc particle velocities with respect to the de- 
tector induces a time dependence in the differential detection rate, i.e. S{Er, t) = 
So{Efi) + Sm{Ef{) cos[uj(t — to)], where u = 27r/365 days and to = 153 days (June 2"^^). 
So{Er) is the average (unmodulated) differential rate defined in Eq.([T9|)and Sm{Eji) 
is the modulation amplitude of the rate. The relative importance of SmiEn) with 
respect to So{Eji) for a given detector, depends both on the mass of the dark matter 
particle and on the value of the recoil energy where the effect is looked at. Typical 
values of Sm{Eji)/ So{Eji) for a Nal detector range from a few percent up to ~ 15%, 
for neutralino masses of the order of 20-80 GeV and recoil energies below 8-10 KeV. 

The search for the annual modulation effect is currently undertaken by the DAMA/ 
Nal Collaboration^. The analysis of their data over two years of data-taking pro- 
vides the indication of a modulated signal 0. This result can be shown as the closed 
contour in the ^cr^caiar'"'^'' pl^^i^ displayed in Fig. 2. The region inside the 

contour is compatible with a modulation signal at 2-cr level. The contour takes into 
account also the uncertainties in astrophysical velocities, as discussed in 0. Fig. 2 
shows that there exist susy configurations which would be able to explain such an 
effect. In the papers of Refs.BS it has been proved that some of these configurations 
are explorable at accelerators and/or by WIMP indirect experiments and that the 
relevant relic neutralinos might behave as major components of cold dark matter. 
For an analysis of these configurations in a SUGRA scheme, see also Ref.E^. 



6.4- Indirect detection: neutrino fluxes from Earth and Sun 

Other ways of detecting dark matter neutralinos rely on the possibility to detect 
the products of neutralino annihilations. One perspective is to observe a neutrino 
signal coming from celestial bodies, namely Earth and Sun, where the neutralinos may 
have been captured and accumulated during the lifetime of the macroscopic body. The 
sequence of the physical processes which could produce these signals are: i) capture 
of relic neutralinos by the macroscopic bodies; ii) subsequent accumulation of these 
particles in a region around the centre of these celestial objects; iii) pair annihilation of 
the accumulated neutralinos which would generate, by decay of the particles produced 
in the various annihilation final states, an output of high-energy neutrinos. Since the 
process of annihilation takes place inside a medium (i.e., the interior of the Earth 
or the Sun), the annihilation process is perturbed as compared to the annihilation 
in vacuum. This effect can be effectively taken into account0 by neglecting the 
contributions of the light quarks directly produced in the annihilation process or 
in the hadronization of heavy quarks and gluons, because these light particles stop 
inside the medium (Sun or Earth) before their decay. For the case of the Sun, also the 
energy loss of the heavy hadrons in the solar medium and the energy loss of neutrino 
themselves have to be considered^. 
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The differential neutrino ffux is then calculated as 



dE, Aird^f-^ dE, ^ ' 

where Ta denotes the annihilation rate, d is the distance of the detector from the 
source (i.e. the center of the Earth or the Sun), F denotes the neutralino pair an- 
nihilation final states, B^^ denotes the branching ratios into heavy quarks, r lepton 
and gluons in the channel F. The spectra dNf^/dEi, are the differential distributions 
of the neutrinos generated by the r and by hadronization of quarks and gluons and 
the subsequent semileptonic decays of the produced hadrons. A detailed calculation 
of these spectra is usually performed by means of a Montecarlo simulationE3. The 
spectra due to heavier final states, i.e. Higgs bosons, gauge bosons and t quark, can 
be computed analytically by following their decay chain down to the production of a 
b quark, c ouark or a tau lepton, where the result of the Montecarlo simulation can 
be appliedffl. 

The neutrino flux may be detected in a neutrino telescope by measuring the muons 
which are produced by and z/^ interactions in the rock around the detector and 
then traverse it upwardly. Therefore, the signal consists of a flux of up-going muons, 
which is computed as 

r da^E^lJdN, 

where X is the muon range in the rock, da^^E^,, E'^)/dE'^ is the charged current 
cross-section for the production of a muon of energy E'^ from a neutrino of energy E^, 
and Psurv{Efj,, E'^; X) is the survival probability that a muon of initial energy E'^ will 
have a final energy E'^ after propagating along a path-length X inside the rock which 
surrounds the detector. The function Psurv(-£'/i, E'^] ^) therefore takes into account the 
energy losses of muons in the rock. In Eq. (piD , E^^ is the detector threshold energy, 
which for current neutrino telescopes like MACRO and Baksan is of about 1-2 GeV, 
and this is quite suitable for neutralino indirect detection, especially for neutralinos 
lighter than about 100 GeV. Large-area neutrino telescopes with higher threshold 
energies (above a few tens of GeV) are more suitable for heavier relic particles. 

Experimentally, one searches for a statistically significant excess of up-going 
muons above the muon flux originated from atmospheric neutrinos. The different 
angular behaviour of the signal with respect to the atmospheric neutrino background, 
which has a rather flat distribution as a function of the zenith angle, is the handle 
for the signal-to-background discrimination. Clearly, the flux from the Sun can be 
pointed at directly towards the direction of the Sun. In the case of the flux from the 
Earth, the process of accumulation of neutralinos induces a rather peaked distribution 
of the neutrino source around the Earth's center. Indeed, the angular distribution is 

G{e) ~ Am^a exp{-2m^a sin^ 6) (22) 
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where 6 is the zenith angle and a = 1.76 GeV^^. This means that for neutrahnos 
(which are heavier than ~ 25 GeV) the signal is produced inside a region whose 
angular extension is less than about 10 degrees. 

The experimental searches at neutrino telescopes have found no muon excess so 
far and therefore upper limits on the muon flux have been set. The solid line in 
Fig. 3a is the current most stringent experimental upper limit from the MACRO 
CoUaborationSii for the neutrino flux from the center of the Earth. Fig. 3b refers 
to the flux from the Sun. 

Again, superimposed to the experimental limits, we show in Fig. 3a (Earth) and 
Fig. 3b (Sun) the susy scatter plot for the up-going muon signal ^^^^^^ and 
The susy configurations have been varied in the ranges quoted in Sect. 6.1. The 
astrophysical parameters are pi = 0.3 GeV cm~^, vq = 220 Km s~^, Wesc = 550 Km 
s~^ and the cosmological bound Q^h"^ < 0.7 has been applied. 

6.5. Indirect detection: antimatter and gamma rays 

The annihilation process of dark matter neutralinos may take place also directly 
in the galactic halo. In this case, many different signals other than neutrinos are 
possible. These signals, at variance with the signals previously discussed, which 
depend only on local galactic properties, depend directly on the matter distribution 
of neutralinos over the whole Galaxy. Moreover, the propagation inside the Galaxy 
of the particles which constitute the signal is perturbed by the Galaxy itself, like, for 
instance, interactions with the interstellar medium or, in the case of charged particles, 
diffusion in random magnetic fields. 

One of the most interesting possibility is the production of antimatter from neu- 
tralino annihilation in the halo. The fluxes have been calculated for production of 
antiprotons, antideuteron and positrons. 

Antiproton^^^ and Antideuteron^ can be produced by the decay and hadroniza- 
tion of the final state particles of the annihilation process. The differential rate per 
unit volume and unit time for the production of p's from neutralino pair annihilation 
is defined as 

\ '"'X J F,h "-^P 

where < aann"^ > denotes the average over the galactic velocity distribution function 
of neutralino pair annihilation cross section Uann multiplied by the relative velocity 
V of the annihilating particles, Tp is the antiproton kinetic energy and -B^/^ is the 
branching ratio for the production of quarks or gluons h due to the decay of the 
particles produced by neutralino annihilation into the final state F. Finally, dNp/dTp 
is the differential energy distribution of the antiprotons generated by hadronization of 
quarks and gluons. Notice that the rate depends on the square of the mass distribution 
function of neutralinos in the galactic halo '^^^ '^^^^ production of D is clearly 

analogous to Eq. (^5]) . 
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After being produced, antimatter propagates inside the Galaxy and it experiences 
both diffusion in the galactic magnetic field and energy losses, due to ionization, scat- 
tering, collision and others. The propagation of antiprotons in the galactic medium is 
properly calculated in a diffusion model where the Galaxy is described as composed 
of two zones, one for the disk and the other for the halo. The diffusion equation which 
governs the behaviour of antiprotons is 

V-{kV^p) - 2M(z)rpVp + 2M(^)gr"W - ^ {K^)V^p} = , (24) 

where ipp denotes the antiproton density, K is the diffusion coefficient, h is the height 
of the galactic disk and Fp is the collision probability with the interstellar medium 
and b{E) describes energy losses. An analogous equation holds for the antideuterons. 

Solution of the diffusion equations gives the antiproton (or antideuterons) flux at 
the heliosphere boundaries (interstellar flux) as 

M&.T,) =<a_.> ^ f£^)%f(0,Tp)-E4I^^' (25) 

where '0p^(©, Tp) is obtained by solving Eg . (12^) . 

Antimatter subsequently enters the heliosphere where it propagates against the 
solar wind before arriving at the Earth. The effect induced by the solar wind (so- 
lar modulation) is quite important at low kinetic energies, and introduces a time 
dependence into the calculation, since it is correlated with the 11 year solar cycle. 

Both antiprotons and antideuterons can be produced also by the interaction of 
primary cosmic rays on the interstellar medium. This secondary antimatter fluxes 
constitute a background for the susy signal. The background fluxes have a different 
energy behaviour as compared to the the ones of susy origin. In particular, the low 
energy tail of the energy spectrum is the most interesting place to look at, since the 
signals have a somewhat flat behaviour, while the secondary fluxes are depressed by 
kinematical reasons B This effect is stronger for antideuterons than for antipro- 
tons. Therefore, the antideuteron flux presents some advantages of discrimination 
over background with respect to the antiproton flux, even if its smallness makes the 
detection harder to be achieved. 

Fig. 4 show the antiproton flux0 $p vs. the neutralino mass for the susy conflgu- 
rations introduced in Sect. 6.1. The flux is calculated for ap kinetic energy Tp = 0.24 
GeV, to conform to the flrst BESS energy bin, and for a solar modulation phase close 
to the BESS data-taking periods. The horizontal line is the present upper limit de- 
rived from the BESS 95 and BESS 97 dataS. We notice that antiproton measurement 
are already able to exclude some susy conflgurations which would imply a too large 
p flux at low kinetic energies. 

Fig. 5 shows the antideuteron fluxli vs. the neutralino mass for the same 
susy models. The flux is calculated for a D kinetic energy Tg, = 0.24 GeV/nucleon, 
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at the maximum of solar activity. The horizontal line is the sensitivity level which 
could be achieved by the AMS detector during the flight on space station. 

Positron^ are produced again from the decay chain of the neutralino annihi- 
lation products. It is also possible to produce directly a pair of monocromatic 
electron-positron pair. The branching ratio for this process is usually very small, 
but some susy models can present a production rate strong enough to be at the 
level of the detector sensitivityil. The calculation of the positron flux is analo- 
gous to the one for antiprotons or antideuterons, with the inclusion of additional 
energy loss mechanisms: inverse Compton scattering on starlight and cosmic mi- 
crowave background, and synchrotron radiation emission in the galactic magnetic 
field. The IS positron flux is then affected by the solar wind before coming to the 
Earth where it can be detected. Fig. 6, taken from Ref.E3, shows a scatter plot 
of the positron flux $e+ integrated in the energy range 8.9 - 14.8 GeV, to con- 
form to one of the HEAT data bins. The susy parameters have been varied in the 
ranges: GeV < {M^l < 5000 GeV, GeV < < 5000 GeV, GeV < tua < 
10000 GeV, 100 GeV < mo < 30000 GeV, -3 < A < +3, 1 < tan/3 < 60. Only 
configuration with 0.025 < f^^/i^ < 1 are shown on the plot. The horizontal dashed 
line denotes the value of the background positrons of secondary origin calculated in 
Ref.0. The horizontal band represents the HEAT 94 dataE3. 

The last possibilitv weconsidfir here is the production of diffuse gamma rays or a 
gamma line (see Refs.E^KElHEaJHJ and references quoted therein). Diffuse photons are 
produced mainly through the decay of neutral pions originated from the hadronization 
of the neutralino annihilation products. A monochromatic gamma line, instead, is 
produced through the loop-processes XX — *■ 77 XX — *■ ^7- this case, 
the gamma line would constitute a particularly nice signature, since it is practically 
background free. 

In both cases, the fluxes are usually rather low and, in order to have fluxes at the 
level of the detector sensitivities, some matter over-density is needed, like for example 
a singular dark matter halcS'S or a clumpy matter distribution^. Fig. 7, taken from 
Ref.l^, shows a scatter plot of the gamma ray flux $cont,7 integrated above = 1 
GeV. The susy parameters have been varied in the same ranges defined above for the 
positron flux. Only configurations with 0.025 < Q^h'^ < 0.5 are shown on the plot. 
The horizontal line denotes the integrated gamma ray flux measured by egretH. 
Fig. 8 from Ref. ii shows the perspectives of a number of Air Cherenkov Telescopes 
for the measurements of the expected 7-ray line. 



7. Conclusions 



As we have seen above, the question of dark matter in the Universe presents a 
large number of intriguing facets of relevance for cosmology, astrophysics and particle 
physics. It represents a field in a very fast expansion, because of an impressive 
development in experimental activities as well as in theoretical investigations. Let 
us mention just some of the most promising avenues: a) Further observations and 
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analyses of high-z Supernovae, of the Cosmic Microwave Background Radiation, and 
of the time evolution of the number density of clusters are expected to provide more 
conclusive information on Qm and b) New numerical simulations of cosmological 
structures should give a unique information about the (hot/cold) composition of dark 
matter and about crucial details on the dark matter density profile in single galaxies, 
c) Further accumulation of data in WIMP direct detection are expected to play a 
fundamental role in the process of the identification of dark matter particle candidates. 
Present experimental results have already allowed the pinning down of a sector of the 
supersymmetric parameter space, part of which can be explored at accelerator and by 
WIMP indirect searches. It has been proved that the relevant relic neutralinos might 
behave as major components of cold dark matter. No doubt that the connection of 
particle physics with the dark matter problem in the Universe is one of the most 
exciting and far-reaching field in astroparticle physics. 
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FIG. 1. Neutralino relic abundance Q^h'^ as a function of the neutralino mass m^, 
evaluated in the MSSM. Only configurations which provide a rehc abundance not in 
contrast with the age of the Universe (i.e. Q^h'^ < 0.7) are shown. 
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FIG. 2. Neutralino-nucleon cross section CTg^j^^'^ times the fractional amount of 
neutralino dark matter C,, as a function of the neutrahno mass m^. The figure refers 
to the value pi = 0.3 GeV cm~^ for the total local dark matter density. The solid 
line is the present upper limit at 90% C.L.il. The closed contour delimits the region 
compatible with the annual modulation effect^lH. The scatter plot shows the quantity 
C,cricailr°^^ evaluated in the MSSM. Different neutralino compositions are shown with 
different symbols: crosses for gauginos, open circles for higgsinos and dots for mixed 
neutralinos. 
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FIG. 3a. Flux of up-going muons from neutralino annihilation in the Earth, 

plotted as a function of m^. The solid line denotes the present upper limitS. Different 
neutralino compositions are shown with different symbols: crosses for gauginos, open 
circles for higgsinos and dots for mixed neutralinos. 
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FIG. 3b. Flux of up-going muons ^^'^ from neutralino annihilation in the Sun, 
plotted as a function of m^. The solid line denotes the present upper limitS. Different 
neutralino compositions are shown with different symbols: crosses for gauginos, open 
circles for higgsinos and dots for mixed neutralinos. 
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FIG. 4. Antiproton flux $p as a function of rriy. for a p kinetic energy Tp = 0.24 GeV 
and for the solar minimum. The horizontal line is the present upper limit derived 
from the BESS 95 and BESS 97 date""" 
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FIG. 6. Positron flux integrated in the energy range 8.9 - 14.8 GeV vs. m^, for 
configuration with 0.025 < ^^^^^ — "^^^ horizontal dashed line denotes the value 
of the background positrons of secondary originEll. The horizontal band represents 
the HEAT 94 data. 
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FIG. 7. Gamma ray flux $cont,7 integrated above = 1 GeV vs. m^, for configu- 
rations with 0.025 < ^x^"^ — '^^^ horizontal line denotes the integrated gamma 
ray flux measured by EGREtH. 
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FIG. 8. Gamma ray flux inside a 10 ^ sr angular cone which contains the galactic 
center, for the density matter distribution of Ref.§. The flux refers to the contribu- 
tions of the 27 and 7Z annihilation lines for heavy neutralinos. The solid lines denote 
the estimated sensitivities for different atmospheric Cherenkov detectors. 
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